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Heavy quarkonium production in the framework of the non-relativistic quantum chromodynamics 
and leading order of the parton Reggeization approach at the Tevatron and LHC is discussed. 
In this note, we compare our predictions for the bottomonium production at the LHC due to 
the color-singlet approximation of the non-relativistic quantum chromodynamics with CMS and 
LHCb data. It is found, that in the production of T(15') states, the color-singlet mechanism is 
dominating, whereas to describe the data for the inclusive T(25) and T(35') production, the color- 
octet contributions should be taken into account. 
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1. Introduction 

The production of heavy quarkonium at hadron colliders provides a useful laboratory for test- 
ing the high-energy limit of quantum chromodynamics (QCD) as well as the interplay of pertur- 
bative and nonperturbative phenomena in QCD. This high-energy processes are dominated by the 
multi-Regge final states, when the contribution of partonic subprocesses involving t-channel parton 
exchanges to the production cross section can become dominant. Then the transverse momenta of 
the incoming partons and their off-shell properties can no longer be neglected, and we deal with 
"Reggeized" ^-channel partons. 

We use gauge-invaiiant approach to introduce off-shell properties of the incoming partons to 
the matrix element of the hard subprocess, which is based on the gauge-invariant effective action 
for QCD in the Regge limit [|l], ^. This effective field theory incorporates besides ordinary gluons 
and quarks, also fields of the Reggeized gluons and quarks. CoiTcsponding Feynman rules where 
obtained in Refs. |^ ||, ^], and they can be used to obtain amplitudes with off-shell incoming 
partons. 



2. Heavy quarkonium production in the LO PRA. 

In the present note, we study the process of heavy quarkonium production, exploiting the 
leading order (LO) of parton Reggeization approach (PRA) and non-relativistic QCD (NRQCD) 
factorization hypothesis [^, ^]. Working at the LO in a, and the relative velocity of heavy quarks 
V, we consider the following partonic subprocesses, which describe heavy quarkonium production 
at high energy: 

R{q,)+R{q2) ^ ^^['P\'\'sf\'s^^\'pf\p), (2.1) 

R{qi)+R{q2) ^ J^['sWp) +g{p'), (2.2) 

where R is the Reggeized gluon and g is the Yang-Mills gluon, respectively, with four-momenta 
indicated in parentheses, J^[n\ is the physical charmonium state, produced through the intermedi- 
ate Fock state of the qq pair [^^^^Lj ' ] with the spin S, total angular momentum J, orbital angular 
momentum L and in the color-singlet '^^ (CS) or in the color-octet ^^' (CO) quantum states. The 
squai-ed amphtudes \£/{R + R^ qq[^^+^L) ' '])p of the relevant subprocesses in the LO PRA have 
been obtained in the paper [^], and they are shown to be coinciding with our earlier results [|| in 
the ^r -factorization approach [^. 

In the high-energy factorization, the hadronic cross section is presented as convolution of 
off-shell partonic cross-section and unintegrated parton distribution functions (uPDF). The best 
agreement with experiment for the charmonium and bottomonium production, as well as for other 
processes at low pj [1T0|], was obtained with the use of the KMR | [TT| ] procedure, which prescribes 
how to obtain the uPDFs from the conventional integrated ones. For the more details of the for- 



malism, computational formulas and the fitting procedure for the CO NMEs, see the Ref. [12] and 
references therein. 
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Figure 1: Comparison of the LO PRA CS model predictions for the inclusive T(«5') production with the 
CMS jig ] (a), and LHCb [[IJ] (b, c) data. Dashed lines - direct production, dash-dotted lines - feddown 
decays contribution, soUd lines - their sum. The branching -fractions B{Y{nS) -^ fi^jj.^) are included. 

3. Numerical results. 



Our latest numerical results for the prompt charmonium production in the Regge Umit of QCD 



are described in the Ref. []12p. We have shown that it is possible to describe the LHC experimental 
data on the prompt charmonium production at the y^ = 7 TeV in a wide kinematical range (2 < 
Pt < 20 GeV and \y\ < 3.5) with a good accuracy, using the CO NMEs extracted from Tevatron 
data at the \/5 = 1.8 TeV and 1 .96 TeV. The fitted CO NMEs are also shown to be compatible with 



NLO collinear parton model results of Ref. [|13||. 

Here we discuss our results for the inclusive T{nS) transverse-momentum spectra in the LO 
PRA and CS approximation of the NRQCD. As in the Ref. [|l2[], here we used the KMR prescription 



with the MRST-2006 [ 14] set of the PDFs as the collinear input. The values of the CS NMEs were 



used the same as in the Ref. [0], taking into account corrected values of branching fractions for 
the radiative decays of the T{nS) and Xhj{nP) states, accordingly the latest version of Particle Data 
Group [|l5|]. Due to the lack of space, we do not present in this note our new predictions for the old 
Tevatron data which show the qualitative features same with the features indicated below for the 
LHC data. 

In the Figure |l|, the comparison of the LO PRA prediction on the T{nS) production through 
the CS mechanism with LHC experimental data for inclusive T(«5') production in the pp-collisions 
at the \/S = 7 TeV, measured by the CMS [|l^] and LHCb [ |l7| ] collaborations, is presented. The 
factorization and renormalization scales in this calculation are chosen to be /!« = Mf = ^Mt where 
Mj = \/m^ + pj. is the heavy quarkonium transverse mass, and ^ is varied in the interval 1 /2 < 



(^ < 2 to estimate the scale-choosing uncertainty of our calculation, and the result of this variation 
is indicated on the plots by the gray band. 

Figure |I| shows that the CS contributions are not sufficient to describe the data for the pro- 
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duction of the T(25') and T(35) states, whereas the data on the T(15) production in the region of 
Pt < 15 GeV are described quite well. Moreover, description of T(15) production in the rapidity 
interval 3 < |y| < 3.5 becomes better, because of the decrease of the fraction of feeddown-decays 
contribution with rapidity. These features strongly indicate the necessity to include the CO in- 
termediate state contributions to the calculations. As in the case of the chamionium production 



3c(8) 



[g, |12|, |13|], the leading contribution here will be from the CO S\ intermediate state, which px- 
spectrum has the much lower slope, than the slope of the CS contribution. So, as a result of the fit, 
this contribution should affect the inclusive T(15') /77-spectrum only for the pj > 15 GeV, while 
it's contributions to the production of T(25'), T(3S) and f -wave states will be greater. The fit of the 
CO NMEs for bottomonium production will be the subject of our future study. 
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